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excess in terms of dark matter. 
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1. Introduction: the Galactic Center gamma-ray excess 

An excess of gamma rays from the Galactic Center (GC) has been detected in the Fermi-LAT 
data between roughly 0.1 and 10 GeV [1, 2, 3,4, 5]. This excess has a spatial extension smaller than 
10° x 10°, and is spherically symmetric. It was obtained by subtracting to the data known sources 
and a template for the background diffuse emission provided by the Fermi Collaboration. Although 
this background modelling procedure has been debated, the picture that has emerged seems to be 
robust. There is a variety of astrophysical explanations for the excess, but an interpretation in terms 
of dark matter (DM) is nevertheless possible. 

The best fit to the excess quoted in the literature has been obtained for 30 GeV DM annihilating 
into bb (Fig. 1, top left panel), with a cross section of 2 x 10 26 cm 3 s -1 . The data also point to a 
density profile p A 1,2 . A mixture of 90% leptons and 10% b quarks gives a relatively good tit 
(Fig. 1, top right panel), while the tit is bad for a final state containing 100% leptons (Fig. 1, bottom 
panel). Flere the term “leptons” refers to democratic annihilation into leptons, i.e. a combination of 
the e + e~ , p 1 p A r 1 t final states, with 1/3 of the annihilations into each of these channels. 

Consequently, final states of DM annihilation containing only leptons have not been consid¬ 
ered viable when interpreting the GeV excess in terms of DM. However, these conclusions were 
obtained taking into account only the prompt gamma-ray emission, namely the final-state radiation 
(FSR) single-photon emission, and the immediate hadronization and decay of the DM annihilation 
products into photons. Nevertheless, electrons and positrons are also by-products of DM anni¬ 
hilations, and they produce gamma rays via inverse Compton (IC) scattering off photons of the 
interstellar radiation field and bremsstrahlung. It had been argued by the authors of Refs. [6, 7] 
that the contributions of these secondary emissions to the gamma-ray spectrum should not be ne¬ 
glected and should lead to corrections. On top of that, diffusion of electrons and positrons through 
the Galactic magnetic field must be taken into account when modelling these emissions. Here we 
show that these contributions do not just give corrections but they totally change the interpretation 
of the excess in terms of DM. 

2. Diffusion of electrons and positrons from DM 

In order to compute the gamma-ray spectrum from DM-induced electrons and positrons, we 
solve the diffusion-loss equation of cosmic rays. Assuming a steady state, the equation reads [8] 

^ v V+^(VV / ') + f/ = °, (2.1) 

where \\f(x, E) is the cosmic-ray spectrum after propagation, K is the diffusion coefficient: K(E) = 
Kq (£/£()) with Eo = 1 GeV, b tot (E) is the total energy loss rate (IC, synchrotron, bremsstrahlung) 
and q(x,E) is the source term for DM annihilations, proportional to p 2 . To solve the equation, we 
use the semi-analytic method described in Ref. [8]. In this approach, the spectrum of electrons and 
positrons after diffusion is given by 

k* r°° h ii 

xis(x,E) = -— UME,E s ))—(E s )dE s , (2.2) 

b tot {E) Je d E 
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Figure 1: Best fit to the spectrum of the residual extended emission in the 7° x 7° region around the GC, 
for 30 GeV DM annihilating into 100% b quarks (top left panel), and 10 GeV DM annihilating into 90% 
leptons and 10% b quarks (top right panel), and 100% leptons (bottom panel). The best-fit cross section is 
~ 2 x 10 26 cm 3 s'. The data points are taken from [3]. 


where K = (1/2) (<rv) (Po/ ,77 dm) 2 > with (av) the annihilation cross section, p 0 the DM density in 
the Solar neighborhood, and /??dm the DM mass. // is the so-called halo function that contains all 
the information on diffusion, through the diffusion length Ad(£. E$) that depends on the injection 
energy it's and the energy after propagation E. The halo function is convolved with the injection 
spectrum dn/dE. The most difficult step of the resolution is to compute the halo function in the 
context of a cuspy DM profile. For that we used the method relying on Green’s functions. The 
halo function is thus given by the convolution over the diffusion zone (DZ) of the Green’s func¬ 
tion G(x,E',xs,Es) = G{x.xv,.X\)(E.Es)) of the diffusion-loss equation and the square of the DM 
density [8] 

h{X D (E,E S ))= [ dx s G(T,i?;T s ,£s)f—) • (2.3) 

7dz V P© / 

The difficulty with this integral is twofold. First of all, to deal with the steepness of the DM profile, 
we used logarithmic steps. Second, the halo function must boil down to (p/p©) 2 when diffusion 
becomes negligible, that is when Ad —> 0 (i.e. E —> E$). The problem is that in this regime, G 
becomes infinitely peaked, and the sampling of the integral must be carefully chosen in order not 
to miss the peak. The solution is to define different regimes for G depending on the ratio of Ad and 
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the distance to the GC, as described in detail in [9]. Once / has been computed with this dedicated 
treatment of diffusion on small scales, one can compute the IC and bremsstrahlung gamma-ray 
fluxes from electrons and positrons. 

3. Fitting the Fermi-LAT GeV excess with the leptonic channels 

We used fixed values for the quantities describing the interstellar medium (the magnetic field 
and the gas density) that correspond to average values for the losses in the GC region. In Fig. 2 
we show that for democratic annihilation into leptons, the contributions from prompt emission, IC 
and bremsstrahlung are of the same order of magnitude, and they combine to give an excellent fit 
to the excess, with a best-fit cross section of (<7v) = 0.86 x 10” 26 cm 3 s . This is a very important 
result, since it means that 30 GeV DM annihilating into bb is not the only possible annihilation 
set-up, and that DM can annihilate into leptons. This result is discussed in more detail in Ref. [10]. 
We also show in Fig. 3 (left panel) that including the contributions from IC and bremsstrahlung 
does not significantly affect the spectrum for the bb channel, except at low energies. 
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Figure 2: Spectrum of the residual extended emission in the 7° x 7° region around the GC. The data points 
are taken from [3]. The prompt (black dashed), IC (green dashed-dotted) and bremsstrahlung (red dotted) 
emissions from 10 GeV DM democratically annihilating into leptons add up to give a very good fit to the 
data, as shown by the black solid line, with a best-fit cross section of (av) = 0.86 x 10~ 26 cm 3 s _1 . 


Therefore, the effect of including the IC and bremsstrahlung contributions is maximal for 
democratic annihilation into leptons. However, the authors of Ref. [11] used the AMS data to 
set constraints on the leptonic channels, which exclude annihilations into e e and impose that 
the branching ratio into ji ' jl should not exceed 25%. We do not discuss the validity of these 
constraints here, but we show in Fig. 3 (right panel) the best fit to the excess obtained with a final 
state of DM annihilation containing 25% and 75% T' T . With such branching ratios, the fit 

is marginally good, with an effect of the secondary contributions to the gamma-ray spectrum less 
significant than for democratic annihilation. 

Finally, a very important test of the leptonic scenario is the morphology of the emission. In¬ 
deed, due to spatial diffusion, the IC and bremsstrahlung emissions are more extended than the 
prompt component, and their morphology depends on the observed energy. This may allow one 
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Figure 3: Best fits to the Fermi residual with the gamma-ray spectrum from annihilations of 30 GeV DM 
particles into 100% bb (left panel), and 10 GeV DM annihilating into 25% p p and 75% T + T~ (right 
panel). The best-fit cross section is ~ 2 x 10 26 cm 3 s _1 for the left panel, and ~ 1 x 10 26 cm 3 s 1 for the 
right panel. 

to set constraints on the annihilation channel. Shown in Fig. 4 (left panel) is the gamma-ray flux 
plotted against latitude, for the best-fit parameters corresponding to the spectrum of Fig. 2, and 
Ey = 0.1 GeV. Above a few degrees the morphology should be compatible with the one found in 
the literature, and corresponding to prompt emission. However, between £?(0.1) and f)'( I)°, the 
flux profile is shallower than that of the prompt component, which might lead to a tension between 
the flux from the leptonic channels and the morphology from the literature. However, 0.1 GeV is 
below the lowest energy data point (around 0.3-0.4 GeV), and the secondary emissions dominate at 
low energy. It turns out that at 1 GeV, the tension between the total flux from the leptonic channels 
and the morphology of the prompt emission is much weaker, as shown in Fig. 4 (right panel). 
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Figure 4: Gamma-ray flux from DM annihilating exclusively into leptons democratically, as a function of 
latitude b, for gamma-ray energies of 0.1 GeV (left panel) and 1 GeV (right panel). 


Conclusion 

The Fermi-LAT GeV excess is a strong case for DM, and annihilations of 30 GeV particles 


5 









































Fitting the Fermi-LAT GeV excess: on the importance of the propagation of electrons from dark matter 
Thomas Lacroix 


into bb provide the simplest set-up a priori. However, it is very important to take into account all 
relevant emission processes and spatial diffusion. Including all these processes drastically changes 
the interpretation of the excess in terms of DM. Therefore, we showed that bb is not the only 
possible channel and DM can in fact annihilate into leptons. Finally, the morphology below ~ 1° 
at low energies can help to discriminate between the leptonic and bb channels. 
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